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I. INTRODUCTION 

When a gas of massive bosons are cooled to a tempera- 
ture very close to absolute zero in an external potential, a 
large fraction of the atoms collapse into the lowest quan- 
tum state of the external potential forming a condensate 
known as a Bose Einstein condensate (BEG) [1-4]. Bose 
Einstein condensation is an exotic quantum phenomenon 
observed in dilute atomic gases and has made a huge 
turnaround in the fields of atom optics and condensed 
matter physics. The recent experimental realization of 
BEGs in rubidium [5] has really kickstarted the upsurge 
in this area of research leading to flurry of activities in 
this direction while the observation of dark [6] and bright 
solitons [7], periodic waves [8], vortices and necklaces [9] 
has given an impetus to the investigation of this singular 
state of matter. 

The dynamics of BEGs is in general governed by an 
inhomogeneous nonlinear Schrodinger (NLS) equation 
called the Gross-Pitaevskii (GP) equation and the be- 
haviour of the condensates depends on the scattering 
length (binary interatomic interaction) and the trapping 
potential. Eventhough the GP equation is in general non- 
integrable, it has been recently investigated for specific 
choices of scattering lengths and trapping potentials [10- 
12]. 

It is known that a BEG comprises of coherent matter 
waves analogous to coherent laser pulses and the wave 
nature of each atom is precisely in phase with that of 
every other atom. In other words, the atoms occupy the 
same volume of space, move at identical speeds; scatter 
light of the same color and so on. Hence, it looked fun- 
damentally impossible to distinguish them by any mea- 
surement. This quantum degeneracy arising out of high 
degree of coherence has been recently exploited in an ex- 
periment by Andrews et al. and Ketterle's group [13] 
showing that when two seperate coluds of BEGs over- 
lap under free ballistic expansion, the result is a fringe 
pattern of alternating constructive and destructive inter- 
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ference just as it occurs with two interesting laser radia- 
tions. In other words, when BEGs were made to collide 
upon release from the trap, de Broglie wave interference 
pattern containing stripes of high and low density were 
clearly observed. These experiments which underlined 
the high degree of spatial coherence of BEGs led to the 
creation of atom laser [14]. Gan one observe the same 
matter wave interference pattern by allowing the bright 
solitons which are condensates themselves to collide in a 
trap? Motivated by this consideration, we investigate the 
collisional dynamics of condensates in a time dependent 
trapping potential. 



II. GROSS-PITAEVSKII EQUATION AND 
LAX-PAIR 



At the mean field level, the time evolution of macro- 
scopic wavefunction of BEGs is governed by the Gross- 
Pitaevskii (GP) equation. 
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where g = ATTh?as{t)/m, V = Vq + Vi, Vo{x,y) = 
mu}\{x'^ Vi = TOWo(0^^/2- In the above equation, 

Vo and Vi represent atoms in a cylindrical trap and time 
dependent trap along z-direction respectively. The time 
dependent trap could be either confining or expulsive 
while the time dependent atomic scattering length as{t) 
can be either attractive (og < 0)or repulsive (os > 0). 
As a result, the condensates confront with both time de- 
pendent scattering length and time dependent trapping 
potential. 

Considering BEGs as an assembly of weakly interacting 
atomic gases, eq. (1) can be effectively reduced to a one 
dimensional GP equation taking the following form (in 
dimensionless units). 
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where ^{t) = —2as{t)/aB, M(t) = ujQit)/^^, as is the 
Bohr radius, M(t) describes time dependent harmonic 
trap which can be either attractive {M{t) > 0) or ex- 
pulsive {M{t) < 0). Now, to generate bright solitons 
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of eq. (2) for both regular and expulsive potentials, we 
introduce the following modified lens transformation 



tPiz, t) = y^A(t)Q{z, t)exp(i$(z, t)) (3) 
where the phase has the following simple quadratic form 



^z,t) = --c{t)z^. 
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Substituting the modified lens transformation given by 
eq. (3) in eq. (2), we obtain the modified NLS equation 

iQt + ^Qzz-ic{t)zQ,-ic{t)Q+'y{t)A{t)\Q\^Q = 0, (5) 
with 



and 
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<t) = -^lnA(t). 
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Equation (5) admits the following linear eigen value 
problem 
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In the above linear eigenvalue problem, the spectral 
parameter which is complex is nonisospectral obeying the 
following equation 

C{t) = cmt) (10) 

with j{t) = 1/A(t). It is obvious that the compatability 
condition {4iz)t = {'t't)z generates eq. (5). 

Substituting eq. (7) with 7(f) = in eq. (6), we 

get 

7"(i)7(i) - 27'(t)' - M{t)^{tf = (11) 

Thus, the solvability of the GP eqn.(2) depends on the 
suitable choices of scattering length 7(f) and the trap 
frequency M(t) consitent with eq. (11). Recently eq. 

(2) has been investigated and bright solitons have been 
generated employing gauge transformation approach [11]. 



III. BRIGHT SOLITON INTERACTION AND 
MATTER WAVE INTERFERENCE 

To investigate the collisional dynamics of condensates 

in the presence of a time dependent trap, we now consider 
the two bright soliton of the following form, 
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FIG. 1: The variation of scattering length 7[f]=0.2sec[-s/2t] 
corresponding to case (i). 



where 



= {-2p2[{a2-a,f -{131-131)] 
- 4i/3i/32(a2-ai)}e('*i+'«^) 
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To generate the interference pattern, we now allow the 
bright solitons to collide with each other in the presence 
of trap for suitable choices of scattering length ^{t) and 
trap frequency M(t) (or c(t)) consistent with eq. (11). 

Case (i); When c(t)=\/2tan\/2i, the variation of scat- 
tering length 7(i)=0.2sec-\/2t is shown in fig.l. Ac- 
cordingly, the trap frequency M(t) becomes a constant 
(M(t)=2). Under this condition, the collisional dynam- 
ics of bright solitons in the harmonic trap is shown in 
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FIG. 2: Interaction of two bright solitons forming interference 
pattern for the choice c[t]=\/2tan[y2t], 7[f]= 0.2sec[-\/2t], 
aio=.01, Q2o=0.8, /?io=0.06, /32o=0.012, 7o=0.2, <j)^ =0.01, 
(j>2=0.1, 5i=0.1, 52=0.01 and (b) the corresponding contour 
plot depicting matter wave density evolution 
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FIG. 3: The phase change evolutions of fig. (2a) 



fig. 2a and the corresponding density evolution in fig. 2b. 
The density evolution consists of alternating bright and 
dark fringes of high and low density respectively while 
the phase difference between the condensates shown in 
fig. 3 continuously changes with time. 

Case(ii): When c(t)=A/2cot\/2i, the scattering length 
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FIG. 4: The variation of scattering length 7[t]=0.2sin['y2f] 
corresponding to case (ii). 

(a) 
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FIG. 5: Two soliton interaction forming interference pattern 
for the choice c[t] = ^cot[V2f], ■y[t]^0.2 sm[y/2t] with the 
other parameters as in fig (la) and (b) the corresponding 
contour plot depicting matter wave density evolution. 



varies sinusoidally with time as in the case of Fesh- 
bach resonance (shown in fig. 4) while the trap becomes 
M{t)=-2cot^V2t - 2cose(?y/U. The evolution of the 
condensates in the above time dependent trap and the 
corresponding density evolution are now shown in figs 5a 
and 5b. From the interference pattern, one observes that 
the matter wave density is maximum at the origin and 
it decreases gradually on either sides. Again, the phase 
change between the condensates oscillates with time as 
shown in fig. 6. 

Thus, our results reinforce the fact that the matter 
waves originating from the condensates (bright solitons) 
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FIG. 6; The phase change evolution of fig. 5a 

do interfere and produce a fringe pattern analogous to 
the coherent laser beams and the interference pattern is 
a clear signature of the long range spatial coherence of 
the condensates. It should be mentioned that the inter- 
ference patterns were obtained earlier by Andrews et al. 
[13] under the condition of free ballistic expansion (atoms 
are released from the trap) while we selectively tune the 
frequency of the trap M(t) in accordance with the scat- 
tering length 7(t) (consistent with eq.(ll)) to observe 
matter wave interference in the collisional dynamics of 
bright solitons. It should be mentioned that the optical 
traps have opened up the possibility of realizing different 



types of temporal variation of the trap frequency M(t) 
while the scattering length ^{t) can be controlled both by 
Feshbach resonance as well as through the trap frequency 
M(t). The phase change evolution which gives a measure 
of the coherence of the condensates shown in figs. 3 and 6 
shows that it oscillates with time and this is reminiscent 
of Josephsen effect wherein the phase difference between 
the superconducting wave functions oscillates with time. 



IV. CONCLUSION 

In conclusion, we have shown that it is possible to ob- 
serve matter wave interference pattern by studying the 
collisional dynamics of condensates (bright solitons) for 
suitable choices of scattering length ^{t) and trap fre- 
quency M(t) and this could happen without free ballistic 
expansion. We also observe that the phase difference be- 
tween the condensates oscillates with time clearly show- 
ing the manifestation of Josephson effect. It would be 
interesting to interpret the matter wave destructive in- 
terference and whether this would mean that atoms plus 
atoms add up to vaccuum. 
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